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The equilibrium structure, energy bands, phonon dispersions, and s- and d- 
channel electron-phonon interactions (EPIs) are calculated for the infinite- 
layer superconductor CaCuC>2 doped with 0.24 holes per CuOi- The LDA 
and the linear-response full-potential LMTO method were used. In the equi- 
librium structure, oxygen is found to buckle slightly out of the plane and, as 
a result, the characters of the energy bands near ep are found to be similar 
to those of other optimally doped HTSCs. For the EPI we find \ s ~0.4, 
in accord with previous LDA calculations for YBaiCu-zO-j . This supports the 
common belief that the EPI mechanism alone is insufficient to explain HTSC. 
X x 2_ y 2 is found to be positive and nearly as large as X s . This is surprising 
and indicates that the EPI could enhance some other d-wave pairing mech- 
anism. Like in YBa2Cu2,Oj, the buckling modes contribute significantly to 
the EPI, although these contributions are proportional to the static buckling 
and would vanish for flat planes. These numerical results can be understood 
from a generic tight-binding model originally derived from the LDA bands of 
YBaaCu^Oj. In the future, the role of anharmonicity of the buckling-modes 
and the influence of the spin-fluctuations should be investigated. 
PACS numbers: 74.72. Jt, 74.25.Kc, 63.20.Kr 
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1. INTRODUCTION 

The mechanism of high-temperature superconductivity in hole-doped 
Cu02 materials remains a subject of vivid debate.HI The Coulomb repulsion 
between the electrons is poorly screened, at least for low doping levels, and it 
may be responsible for the pairing.0 This is supported by experimental evi- 
dence that the symmetry of the paired state is d x 2_ y 2 with lobes in the direc- 
tion of the CuO-bond. Nevertheless, the results of Hubbard- and i-J-model 
calculations using realistic parameters lead to a suspicion that something 
more than the Coulomb repulsion is needed.! A large amount of experi- 
mental data such as superconductivity-induced phonon renormalizations,! 
a large isotope effect away from optimal dopingja and phonon-related fea- 
tures in the tunnelling spectrafca show effects of the electron-phonon interac- 
tion (EPI). It is therefore of interest to investigate whether the EPI, which 
gives a pair-interaction, V (k, k') oc — \g (k, k')| , which is always negative 
(attractive), could support d x 2_ y 2 -wave pairing with a gap-anisotropy like 
A (k) oc cos (ak x ) — cos {ak y ). For this to occur, |g (k, k')| 2 must be large 
when A (k) and A (k') have the same sign, and small when they have oppo- 
site signs. Such an investigation is the main topic of the present paper. 

We first report ab initio local density-functional (LDA) calculations of 
the crystal, electronic, and phononic structure, and of the s- and (i-channel 
EPIs for the infinite-layer compound CaCu02- Unfortunately, no single- 
crystal measurements exist for an infinitedayer HTSC, and the ceramic sam- 
ples for which a T c of 110K was reported^ may be phase impure.13 However, 
the infinite-layer structure is uniquely simple for theoretical studies and high- 
temperature superconductivity is generally believed to take place in Cu02 
planes. As we shall see, also the LDA electronic structure of this mate- 
rial, properly relaxed, is very similar to that of the stoichiometrically doped, 
well-characterized HTSC, YBa2Cus07. Our calculations are for 24 per cent 
hole-doped CaCu02 with the hole charge being neutralized by a homoge- 
neous, negative background charge. We have chosen a doping of 0.24 holes 
per Cu02 because this places the Fermi level at the uppermost saddle-point 
in the LDA band structure, and because it agrees with the experimental 
(average) doping level. 

In the second part of the paper we present analytical low-energy bands, 
originally derived from LDA calculations for YBa2Cu307,Ei and derive the 
EPI for the buckling-mode. This tight-binding model has previously! been 
used to explain the in-plane dispersion of the inter-plane hopping integral 
t± (k||^, which is a crucial ingredient of the inter-plane tunnelling model 

of HTSC0When augmented with a mean-field treatment of the Coulomb 
repulsion U between two holes in the same Cu d x 2_ v 2 orbital, the model has 
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Fig. 1. Band structure of doped CaCuC>2 with fiat planes. 

correctly predicted the value of the inter-plane exchange coupling J± .0 In 
the present paper we consider the EPI in optimally and overdoped HTSCs 
and shall assume that U is screened completely. We shall derive an analytical 
expression for the coupling g (k, k + q) of the electrons to the buckling mode 
for a single plane and show that, for this mode, A for d x 2_ y 2-wacve pairing is 
positive definite. The analytical tight-binding model should also be useful 
for further studies. 



2. LDA CALCULATIONS FOR CaCu0 2 

2.1. Energy bands and equilibrium structure. 

In this section we discuss the electronic energy bands and the equi- 
librium structure calculated with our full-potential linear-muffin-tin-orbital 
(LMTO) method. The bands presented in Fig. 1 agree in detail with those 
of previous calculations, 111 and they hardly change upon hole doping by 10 
per cent, except for a downwards shift of the Fermi energy to the upper- 
most van Hove singularity (at R) as shown in the figure. The Fermi surface 
(FS) for the doped compound (see Fig. 4a) has a square (lO)-oriented cross- 
section with the rXM-plane [k z = 0) and a (ll)-oriented cross section with 
the ZRA-plane 

Although the band structure in Fig. 1 has antibonding in-plane pda- 
character (O x -Cu x 2_ y 2-Oj / ), like for all Cu02-HTSCs, at and above ep, there 
are strong atypical features just below €f '■ First of all, the uppermost band 
at X is an antibonding out-of-plane pdir band [the O z -Cu xz band at (^,0, 0) 
and the O z -Cu yz band at (0, f ,0)], whereas at R [(f ,0, f) and (0, \ , f)] the 
uppermost band is the usual pda-band. This is due to a relatively strong 
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dispersion, 2t\ cos ck z , of the pdir band in the z-direction, a dispersion 

which is opposite to that, — 2tj_ ^k||^ cosck z , of the pda band. Secondly, 
near A the top of the antibonding in-plane pdir band (O y — Cu xy — O x )nearly 
reaches ef- 

Our total-energy calculations however predict that this infinite-layer 
structure is unstable with respect to the q=0 oxygen out-of-phase buckling 
mode (B2 U ). We have performed frozen-phonon calculations for both un- 
doped and doped compounds and with both experimental and theoretical 
lattice parameters. The results are shown in Fig. 2. We see that a double- 
well potential exists for such buckling motions, when the system is doped. 
The minimum is at 5° if we use the experimental values: a=7.297 a.u. and 
c/a=0.829, and at 7° if we use our theoretically determined values: a=7.112 
a.u. and c/a=0.895. Note that YBa2Cu30y has a 7° static in-phase dimpling 
of the Cu-0 planes, both experimentally and theoretically. 

The buckling reduces the in-plane t ZjXZ =t z>yz and ty )Xy =t XiXy hopping 
integrals and thereby causes the tops of the out-of-plane and in-plane pdir 
bands to move well below €f- For the theoretically stable structure, our 
bands in Fig. 3 and Fermi-surface cross sections in Fig. 4b are simi- 
lar to those of other cuprates. Comparing in detail with the calculated 
FS of YBa2Cu307,i where the dimple causes the saddlepoint at the Fermi 
level (that of the odd plane band) to be bifurcated by Ak to the positions 
± AA;,0,0), we see that the saddlepoint at the Fermi level in CaCu02 
is not bifurcated away from R (^,0, — ) . The saddlepoints well below 
(that of the even plane band in YBa2Cu307 and that for k z =0 in CaCuC^) 
are however bifurcated away from X in both materials. (The difference be- 
tween the saddlepoints at k z =^ and k z =0 in CaCu02 is due to the opposite 
fcz-dispersions of the pda and pdir bands.) 

By comparison of the FS cross-sections in Figs 4a and b, we see that 
the static buckle squares up the (lO)-oriented cross sections and thus cre- 
ates large, flat (100)- and (OlO)-oriented pieces on the FS, centered near 
respectively (jg,0,o) and (o,{g,o) . 

The calculated instability suggests that out-of-plane buckling occurs 
generally in such materials when the Fermi level is lowered to the top of the 
pdir bands. Inclusion of the Coulomb correlations, would presumably cause 
such instabilities to occur at lower doping levels than we predict with the 
LDA because, in a mean-field treatment, the energy of the majority-spin Cu 
d x 2_ y 2 orbital is lowered. The striped anti- ferromagnetic low-temperature- 
tetragonaiiAF-LTT) phases recently found by neutron scattering in g-doped 
La2Cu04ll2l and the local buckling instabilities found in 15 per cent doped 
La2Cu04 by EXAFS0 are presumably examples of this. Only when pinned, 
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Fig. 2. (a) Unstable ¥><iu optical mode involving out-of-plane and out— of— 
phase displacements of oxygen atoms, (b) Change in the total energy asso- 
ciated with the B2 U mode. Curve 1 corresponds to the undoped compound 
calculated with the experimental lattice constants a exp and c exp . Curve 2: 
Doped compound with a exp , c exp . Curve 3: Doped compound with calculated 
lattice constants. 



do the AF stripes seem to be decremental to superconductivity. 

2.2. Phonon spectrum 

We now present the phonon dispersions calculated for the doped 
compound. We use newly developed linear-response full-potential LMTO 
method^ whose accuracy was proven on lattice-dynamical, superconduct- 
ing, and transport properties for a large variety of metals.t3 The theoretical 
lattice constants were used in this calculation since the shallowness of the 
B2 U -well calculated with the experimental parameters casts serious doubts 
on the validity of the harmonic approximation. The technical details of this 
calculation can be found elsewhereJIZI 

From the calculated phonon dispersions in Fig. 5, we first notice that 
the structure is stable. Secondly, we notice that a low frequency (1.8 THz) 
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Fig. 4. Fermi surface cross sections of doped CaCu02 with fiat (a) and 
buckled (b) planes. The (lO)-oriented cross section is with the TXM plane 
while (ll)-oriented cross section is with the ZRA -plane. The dashed lines 
show intermediate cross sections. 
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Fig. 5. Calculated phonon spectrum of doped CaCu02 for q z =0. Of the 
branches labelled by solid circles, the lowest is a quadrupolar mode involving 
in-plane rotation of oxygen squares around Cu and z-movement of Ca. The 
two uppermost branches labelled by solid circles involve out-of-plane oxygen 
buckling. The numbers give the s-channel mode EPIs, A Sj! , (q) , and, at the 
top of the figure, the sum over all modes, X s (q) = J2u ^s,u (<?) ■ 

mode exists near the point M. This mode is quadrupolar with the oxygen 
squares rotated around the copper atoms. It can also be viewed as an in- 
plane buckling mode and is strongly coupled to the in-plane pdir bands. The 
mode also involves displacements of the Ca atoms in the z-direction. It is 
interesting to note that for the structure with flat planes, this quadrupolar 
mode is unstable (at M, u =3.7i THz). This is connected with the atypical 
feature near point A of the band structure in Fig. 1. Here again we see that, 
for structural stability of this compound, out-of-plane distortion is required. 

2.3. Linear electron-phonon interaction 

A large amount of theoretical work has been done to estimate the EPI 
in the cuprates using the frozen-phonon approach, i3t3 the rigid- muffin- 
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are numcr- 



tin@ and tight-binding approximations. cJl In addition, there _ 
ous attempts to deduce the strength of the EPI from experiments.^ So 
far, no general agreement is reached, although the LDA frozen-phonon 
calculationst3'll3 for selected, high-symmetry phonons in YBa2Cu30y and 
optimally doped La2Cu04 indicate that \ s ~0.5-1.5 . Compared herewith, 
our linear-response method is more accurate because it allows calculation 
for arbitrary phonon wave vectors q. We have checked this method on EPIs 
for a large number of classical super conductor st3 and the results indicate 
that, granted the LDA, we can obtain reliable estimates of A for the high-T c 
materials. 

Due to the recent experimental evidence for d x 2_ y 2 symmetry of the 
gap,0 we have calculated the strength of the EPI in a general L-channel. We 
use a standard expression for the mode dependent coupling: 



Av,(q) = ^T^T 0) 



i il,v (q ) 

ttN l u,2(q) 
-^-^EM e (k))^(k)|^(k,k + q) 



xY L (k + q)5(e(k + q)) 
where J2k = Ibz d?k/BZV is the average over the Brillouin zone, 

N L ^Y, S(e(k))\Y L (t)\ 2 (2) 

k 

is the electronic 'L-density of states' per spin at the Fermi level (ep = 0) , 
1l,v (q) is the phonon line-width due to EPI in the L-channel, uj u (q) is the 
phonon energy, and we have used Rydberg units. For the electronic energies, 
€j (k) , we have dropped the band index j because the FS of CaCu02 has only 
one sheet. The approximation in the second line of (|l]) is, that the bands are 
linear in the range ep ± uj v (q) . This approximation ought to be scrutinized 
in a future work. Finally, the coupling constant A^ = X^q^L,v(q) is the 
sum over all phonon branches v and the average over the Brillouin zone. 

Fig. 5 shows \ S)l> (q) for the different phonon branches as a function 
of q along high-symmetry directions. The EPI is seen to be particularly 
large for the buckling modes (the two uppermost of the branches marked 
by circles). Remember that for flat planes, these contributions to the EPI 
would vanish. It is furthermore seen that the EPI is strongly enhanced 
along rX, and peaks near TX/3 reaching a maximum of about 2.5 . This 
is partly due to nesting corresponding to sliding of the nearly flat (100)- 
or (OlO)-oriented part of the FS on itself, or on the part parallel to it and 
displaced by (0, g^,0) or (^,0,0), respectively. These sliding-nestings cause 
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Fig. 6. Contributions to the electron-phonon coupling in the d x 2_ y 2- 
channel from different phonon branches. See also caption to Fig. 5. 

( e (k)) 5 (e (k + q)) oc lim w= o x'o (l^) / u to be large near the follow- 
ing flat regions of q-space: q = (q x , 0, q z ) , q = (0, q y , ?*) , q = (q x , ^,q z ) and 
^ = (ia' 1y> 1 Z ) ' w ^h <fa and g y considerably smaller than ^ and consider- 
ably smaller than — . This kind of nesting enhancement is well-known in the 
theory of HTSC.H A second reason for the peaking of A S)I/ (q)near TX/3 is 
a particular wave- vector dependence the EP coupling constants g v (k, k + q) 
for the buckling modes; this we shall return to in the next section where we 
consider the tight-binding model. From Fig. 5 we see that the EPI is small 
in other points of the Brillouin Zone. Averaging over 15 q-points in the 
TXM-plane yields: X s (q z =0) = 0.5, and averaging over the ZRA-plane gives: 
X s (q z =^) = 0.2. These values are respectively upper and lower bounds for 
A s averaged over all q. We conclude that X s ~ 0.4 is too small to account 
for the high-temperature superconductivity in doped CaCu02- 

Finally, we discuss our calculation of the d-channel EPIs. Using eqs. 
(JlJ-j2|) we have analyzed all symmetry-allowed channels of superconductivity 
and came to the conclusion that, apart from the s-channel, relatively large 
positive values of are only obtained for pair-function symmetry d x 2_ y 2. 
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Note that if g were a constant, independent of k and q, would vanish for 
L ^ 0. The X x 2_ y 2 v (q)are presented in Fig. 6 and may be compared with 
\ SjU (q) in Fig. 5. Also for X x 2_ y 2 v (q) there is a strong enhancement near 
rX/3 and, in particular, for the buckling modes. It is also seen that near 
the M-point, \ x 2_ y 2 v (q) vanishes for the buckling modes and is negative 
for all other modes. Averaging over q in the TXM and ZRA planes give 
respectively: \ x 2_ y 2 (q z = 0) = 0.4 and ^ x 2 -y 2 (Qz =] j:) = 0.2. These values 
are positive and only slightly smaller than the corresponding values of \ s . 
The electron-phonon coupling is thus highly anisotropic in 24 per cent hole- 
doped CaCu02, and probably in all HTSCs near optimal doping. 



3. ANALYTICAL MODEL 
3.1. Electronic 8-band Hamiltonian 

In the bottom part of Fig. 7 we specify an electronic 8-band orthogonal 
tight-binding Hamiltonian for a single, dimpled CuC>2 plane, which we feel 
is generic for Cu02-HTSC's and contains the relevant degrees of freedom. 
In the upper Dart, we synthesize the eight plane-bands. This Hamiltonian 
was originallyQ derived from the LDA bands of YBa2Cu307 by integrat- 
ing out the very-high-energy degrees of freedom, neglecting the chains, and 
projecting onto one plane (or, for a bi-layer, to the even or the odd linear- 
combination, or, for an infinite-layer, to a single k z )% We have recently^ 
deduced the parameter values of this Hamiltonian for a number of further 
HTSCs and showed that (apart from Fermi-level position) the only stronglv 
material-dependent parameters are the energy e s of the Cu s _( 32 2_ 1 - ) hybrid,™ 
which depends on the distance of apical oxygen from plane copper, and the 
integral of hopping t z d from the Cu x 2_ y 2 orbital to an O z orbital, which is 
proportional to the buckling-angle 5. For simplicity, we assume tetragonal 
symmetry in Fig. 7 as well as in most ofthe following. For a discussion of 
the inter-plane hopping, we refer to RefB 

Shown in the 2nd row and 1st column of Fig. 7, which we shall refer 
to as Fig. 7(2,1), are the four cr-orbitals looked upon from above the plane, 
|y) = Oj,, \d) = Cu x 2_ J/ 2, \x) = O x , and \s) = Cu s _^ z 2_ l y Note that O y 
and O x refer to orbitals on two different oxygens, 02 at (1, 0, tan 5)a/2 and 
03 at (0, l,tan<5)a/2. For some purposes, we shall also distinguish the two 
oxygens by labels a and b. In Fig. 7(3,1), we see the \d) orbital and two 
of the 7r-orbitals, \z) = O z and \xz) = Cu xz from the edge of the plane. 
Approximate orbital energies (in eV and with respect to the energy of the 
Cu x 2_ y 2 orbital) are given at the relevant points of the band structure shown 
above, in Fig.7(l,l). Hence, e x =e y = e p is 0.9 eV below e^, and e s is 6.5 
eV above. Moreover, e za =e z b = e z is 0.4 eV above e^, and e xz =e yz is 1 eV 
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Fig. 7. The 8-band Hamiltonian for a single CuC>2 plane and synthesis of its 
band structure. The 1st (left) column shows uncoupled pda (O x -Cu x 2_ y 2- 
O y ), Cu s , and pdir (Cu X2 -0 2 and Cu^-O^, stippled) bands. In the 2nd 
(middle) column, the coupling (t sp ) of the Cu s orbital to the pda band has 
been included. In the 3rd (right) column, also the coupling (t z d) between a- 
and 7r-bands has been included. All energies are in eV. 
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below. If we now, as shown in Figs. 7 (2,1) and (3,1) include the pda hopping 
integrals t x d=t y d = t p( f=1.6 eV and the pdir hopping integrals t ZtXZ =t z>yz =0.7 
eV, we obtain the band structure shown in Fig. 7(1,1): The a-orbitals 
give rise to a bonding, a non-bonding, and an anti-bonding O x -Cu I .2_ ?; 2-O y 
band plus a Cu s _(3 2 2_ 1 ) level (full lines), and the 7r-orbitals give rise to two 
decoupled pairs of bonding anti-bonding bands which disperse in either the 
x or y direction (stippled lines). The anti-bonding pda band, which will 
develop into the conduction band, has saddle-points at X (and Y) which are 
well above the top of the 7r-bands, and which are isotropic in the sense, that 
the absolute values of the band masses in the x and y directions are equal. 
This means that the FS at half-filling is a square with corners at X and Y. 
This is shown in Fig. 8(1,1). 

In the 2nd columns of Figs. 7 and 8, the a and ir bands remain decou- 
pled, but we have introduced the Cu s ~O x and Cu s -O y hoppings (t sp =2.3 
eV), as well as the tiny O z -O z hopping (t zz =0.06 eV). The latter is the only 
one reaching beyond nearest neighbors and it merely lifts a degeneracy of 
the anti-bonding 7r-bands at T. We shall neglect it in the following. The 
strong coupling of the remote Cu s _( 32 2_ 1 ) orbital to the pda orbitals has the 
pronounced effect of depressing the conduction band near the saddle-points 
at X and Y, and thereby increasing the mass towards T and decreasing it to- 
wards M (7r, 7r) . This kind of saddle-point we refer to as anisotropic. As seen 
in Fig. 8(1,2), the FS passing through X and Y will therefore correspond to 
a finite hole-doping and will bulge towards T. That the Cu s _( 3z 2_ 1 ) orbital, 
which has the azimuthal quantum number m z =0, may mix with the con- 
duction band at X (and Y), but not along TM follows from the symmetries 
of the corresponding two antibonding pda wave functionsB The net down- 
wards shift of the saddle-point energy is the result of down-pushing by Cu s 
and weak up-pushing by Cu 3x 2_ 1 , whose energy is below the saddle-point. 

As a result of the hybridization with the Cu s _(3 2 2_x) orbital, the saddle- 
points of the anti-bonding cr-band straddle off the top of the appropriate 
7r-band, so that even a weak dimple or buckle of the planes will introduce 
noticeable hybridization between the a and it bands. This is seen in the 3rd 
column of Fig. 7 where we have turned on the weak Cu^^-O^ hoppings 
(*2d=0.24 eV oc sin 5). These couplings become allowed when there is a 
finite angle 8 (« 7° in YBa2Cu30y, YBa2Cu40s, and in calculated CaCuC^) 
between the Cu-0 bond and the plane of the 2D-translations. Since the 
hybridization between the antibonding pda and Cu xz -O z pdn orbitals vanish 
at X (along the XM-line), the a-ir hybridization may make the saddle-point 
bifurcate away from X. This is seen in Fig. 8(1,3). At X (along XM), the 
pda-band does hybridize with the more remote Cu yz -O z pdir band. 
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3.2. Constant-energy contours (CECs) 

The orbitals in the Bloch representation and the 8-band Hamiltonian 
may be found in eq.(2) of Ref.0 and in eq.(l) of RefB As shown in the 
former reference, and given explicitly in its eqs. (12-16), the constant-energy 
contours (CECs), €j (k) = e, have simple analytical expressions. These are 
given in the 3rd and 4th rows of Fig. 8 in order of increasing couplings as 
in Fig. 7. In the present paper we shall neglect the Cu^ and Cu yz orbitals 
and, hence, consider the 6- rather than the 8-band Model (this corresponds 
to using eq.(14) in Ref.i), because then the CECs are just hyperbolas in 
(x, y)-space, where 

x = — (1 — cos ak x ) and y = - (1 — cos ak y ) , (3) 

and they are shown in the 2nd row of Fig. 8. A CEC is specified by 3 numbers 
d, s, and p, which are the values of the scattering functions defined in the 4th 
row of Fig. 8. d (e) specifies the pda interaction, s (e) the spa interaction, 
and p (e) the air interaction. Only d (e) has a non-negligible (increasing) 
energy-dependence in the relevant ±0.1 eV range around ep. For isotropic 
saddlepoints, s — > oo, and for the cases of interest, s ~ 0.5 — 1.5. For a flat 
plane, p = 0. If p = s 2 / (1 + s) 2 , the saddlepoint is extended, that is, the 
dispersion is quadratic towards M, but quartic towards T. If p exceeds this 
value, the saddlepoint bifurcates to the positions given by x = s [pT 1 ! 2 — lj 
and y = 0. It should be noted that, although the CECs near ep can be well 
reproduced by the 6-band model with an energy-independent p, because 
€xz=£yz are several eV below ep and 1.4 eV below e za =e z b, the value of p 
should be adjusted to give the correct CECs. 

The structure of x'o (q><^) obtained_from the 6-band model as functions 
of doping, s, and p is discussed in Ref.E3 

The formalism given above is for tetragonal symmetry and apparently 
only applies to a plane dimpled like in YBa2Cu307, but not to a plane buck- 
led as in Fig. 2. However, since in the 6-band model each Cvorbital couples 
only to the Cu a ,2_ J/ 2-orbital, we can freely define the sign of the O z orbital 
at (1, 0, tan S)a/2, \za) , to be the opposite of that at (0, 1, — tan 5) a/2, \zb) , 
and then the formalism holds also for static buckling. Nevertheless, for the 
sake of preciseness in the wording, we shall assume static dimpling in the 
following. The values of the electronic parameters for various materials may 
be found in RefJ 2 ^ 

In (x, y)-space the average over the Brillouin zone becomes: 




d(e) = (f^fcel s{e) = (££_£K- 



Bifurcated 
Y M 




-d(e) + x + y 

(x-y) 2 

s-\-x-\-y 

+ (2 - x - y)p = 






ic = sin 2 ak x /2 = i (1 — cos ak x ) 
y = sin 2 aky/2 = ^ (1 — cos aA^) 

w = 1 — ic — 2/ = | (cos afcy + cos ak x ) 
v = x — y =| (cos ak y — cos afc^) 
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CAPTION to Fig. 8. The top row of the first three columns shows 
CECs in (k x , k y )-space of the 8-band Hamiltonian corresponding to the three 
columns in Fig. 7. The CECs chosen in the figures are antibonding p<i<7-like 
and are, specifically, the one passing through the saddle-point (e = e sa ddie) 
and two neighboring ones (e = e sa ddic ± Ae). The 4th column illustrates the 
transformation from (k x , k y )-spaee to (x,y)-, or (u,v)-, space. The 2nd row 
shows the CECs in (x, y)-space where they are hyperbolas. The 3rd row 
gives the expressions for the CECs and the (k x ,k y )-(x,y) transformation. 
The 4th row gives the expressions for the scattering functions, and the (u, v)- 
(x, ^transformations. For simplicity, the CECs in the 3rd column are for 
the 6-band Hamiltonian where, of the 7r-orbitals, only the two O^-orbitals 
and not the Cu zx - and Cu„ z -orbitals are included. 



3.3. Coupling of electrons to the buckling mode 

The displacement pattern of the buckling mode is: 

3R 2 <9R 3l z 1 



OQcl ' <9Q q J v / 2MwV / 2 



{exp (-iq • R 02 ) , -exp(-iq-R 03 )} (5) 



where 02 and 03 are the oxygens at respectively (1, 0, tan S)a/2 and 
(0, 1, tan S) a/2 and q is 2-dimensional. 

For q = 0, this is the mode for which the dimple of 02 is increased 
and that of 03 is decreased, or vice versa. For YBa2Cu30y, this is the 330 
cm -1 out-of-phase mode. Had 5 been zero at equilibrium, oxygen would 
have been dimpling upwards in the -Cu-02-Cu- rows and downwards in 
the -Cu-03-Cu- rows, that is what we have called buckling. 

q = (^, 0) : Had 5 been zero at equilibrium, this mode would have been 
buckling in the ^-direction and dimpling in the y-direction. 

q = (-, — ) : Had 5 been zero at equilibrium, both the -Cu-02-Cu- and 
the -Cu-03-Cu-rows would be buckling, and such that the 02-03-02-rows 
in the x=y direction are in phase. This is the displacement pattern of the 
low-temperature orthorhombic (LTO) phase in La2Cu04. 

For this mode, the electron-phonon matrix element is: 

g(k, k' = k + q) = 

— (k| ? 7=2 exp ( "^ q • Ro2) " d^s exp " Ro3 V |k } 



V2M> 
dt zd /dzQ 



R 

c d (k) (cos §fc x | Co (k')| - cos %k y \ c b (k')|) + 
sgn{c a (k') c a (k)} (|c a (k)| cos §*4 - \c (k)| cos %k' y ) c d (k') 
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where, for simplicity, we have only included the phonon-induced modulation 
of the hopping integral t x &. The modulation of e z due to screening is consid- 
ered in RefBs The eigenvector-components, Cd (k) , c a (k) , and q, (k) , for the 
Bloch orbitals |d,k) , \za, k) and \zb, k) are real due to our choice of phases 
for the Bloch functions, and they are most easily obtained by differentiation 
of the CEC-equation with respect to the orbital energies, e^, e za , and e z b. Ac- 
cording to first order perturbation theory, this yields a (k) 2 = de (k) /ded, 
and so on. Simple considerations finally yield the relative signs. As a result, 
we obtain by using the 6-band model: 



S(k,k') 



2t 



zd 



a u a i> 
cos — k x COS —K„ 

2 2 



dt zd /dzp 



cos —k„ cos — k„ 



= c^(l-x)(l-x')-^(l-y)(l-y')y (6) 

where c is the constant in front of the parenthesis on the first line. 

In order to get a feeling for the behavior of the pair-interaction, V oc 
— g 2 , let us first take q = 0, i.e. k' = k. Eq.(|6|) then yields: g = c(y — x) 
which shows that g 2 attains its maximum (c 2 ) for q small and k at those 
points of the FS where the density of states (distance between neighboring 
CECs) is largest and where the superconducting gap is observed to be largest. 
That is, near the saddlepoints. 

Next, let us take k' at the X-point, [x',y ! ] = [1,0] . Then, g = cy/1 — y. 
We thus see that not only does g 2 take its maximal value for k = k' at X, 
but it vanishes for k at X and k' at Y. This is exactly what is needed for 
the electron-phonon interaction to support d x 2_j / 2-pairing. 

Finally, let k and k' to be symmetric around the (1, 1) or (1, —1) line, 
i.e. (k' x ,k'y\ = (ifcy, ±k x ) , in which case q is in the (1,^1) direction. In 

this case, g = V = 0. In particular, this means that V vanishes for q = („j 5) 
where the Coulomb repulsion is maximal. 

Having realized that the buckling mode has the right kind of (k, k')- 
dependence to support d x 2_ y 2-w&ve pairing, let us calculate X x 2_ y 2 for the 
buckling mode and see whether it is really positive. For this purpose we note 
that 

g 2 (k,k') =c 2 x 

(1 - x) (1 - x') + (1 - y) (1 - y') - 2 v /(l-x)(l-y)(l-x')(l-y / ) 
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is factorized so that we can easily calculate \l = X)q Al (q) as a double- 
integral, not over k and q as in ([!]), but over k and k'. Dropping the q- 
dependence of the buckling-mode frequency, we get using (0): 



A. 



4c 2 



1T 2 UJ 



dxdy 



-i 2 



r 1 r 1 s (e) dxdy 

where S (e) = 6 (e (x, y) — £f) ■ Using Y x 2_ y 2 (k) oc x — y we get: 
4c 2 



7T 2 UJ 



dxdy 



n 2 



-y) 



dxdy 



-v) 



Now, the last term of vanishes and we may symmetrize the first 
term as follows: (1 — x) (x — y) — > ^ [(1 — x) (x — y) + (1 — y) (y — x)] = 
— \ [x — y) 2 . As a result: 



1 /"I 



TT^UJ Jo JO 



(x-y) S(e(x,y)) 



dxdy 



(7) 



y/xy(l - x) (1 - y) lo 

which is positive, and easy to evaluate. 

It is obvious from these equations, that if the FS had weight merely at 
X and Y, that is, if were at a non-bifurcated saddlepoint, then 

A s = Xd (-> °o) . 

For €f near an extended saddlepoint, A^ is particularly large. 

The integral (0), as well as the one for the normal density of states, 
N = N s , may be expressed analytically in terms of elliptic integrals if we use 
a band structure with p = 0, i.e. the 4-band model. The results for s = 100 
and s = 0.6 are shown graphically in Fig. 9 as a functions of the doping. 
We see that for doping to the van Hove singularity, = N, as was stated 
above. 



4. CONCLUSIONS 

We conclude that the linear electron-phonon interaction for buckled 
planes may support, but is hardly sufficient to cause, high-temperature su- 
perconductivity based on d x 2_ J/ 2-pairing. The most important mode seems 
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Fig. 9. Density of states and d x 2_ y 2 weighted density of states as denned 
in eq. 2 for the 4-band model, n is the number of electrons counted from 
the half-full band. 

to be the buckling mode, because it modulates the saddlepoints of the en- 
ergy bands where the density of states is high and where the superconducting 
gap is observed to be maximum. Moreover, the buckling mode does not in- 
teract with the electrons for q near (f, f ) where the Coulomb repulsion is 
maximum. 

Considering the facts that the linear EPI for the buckling mode is pro- 
portional to the small static (on the scale of phonon frequencies) buckling, 
and that this mode is calculated to be highly anharmonic, it seems imminent 
to investigate the role of anharmonicity. Similarly, the interplay between the 
electron-phonon interaction and the Coulomb repulsion should be studied. 
For these purposes the analytical 6-band models seem realistic and tractable. 
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